Effective spin model for interband transport in a Wannier-Stark lattice system 
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We show that the interband dynamics in a tilted two-band Bose-Hubbard model can be reduced 
to an analytically accessible spin model in the case ol resonant interband oscillations. This allows 
us to predict the revival time of these oscillations which decay and revive due to inter-particle inter- 
actions. The presented mapping onto the spin model and the so achieved reduction of complexity 
has interesting perspectives for future studies of many-body systems. 
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I. INTRODUCTION 

An amazing control of quantum degrees of freedom 
is nowadays routinely possible with the techniques of 
preparing and handling ultracold matter in the labora- 
tory [lMa]- Backed by a plethora of theoretical propos- 
als (see, e.g., [l-Q), a new direction is the coupling of 
such matter to additional degrees of freedom, such as 
provided by internal states (e.g. [13), by external po- 
tentials (e.g. (HI), by coupling to a bath (e.g. to 
a continuum ( e.g . fl3l flij ]) or even to macroscopic ob- 
jects (e.g. (HjQil). Such hybrid quantum systems are of 
high interest for applications, ranging from fundamental 
physics to metrology. 

A major challenge in studying these systems is to re- 
duce their inherent complexity. This is important for an 
understanding of both the internal dynamics as well as an 
extension to include a coupling to further degrees of free- 
dom. In this paper we focus on the dynamics of atomic 
bosons in a two band Bose-Hubbard model. The prob- 
lem is non-stationary due to an additional Stark force (or 
constant tilt). In particular, we describe the Rabi-like 
oscillations between the two bands, which are well pro- 
nounced in the case of single-particle resonant tunnelling 
between the levels of adjacent lattice wells [l3|, [13 ■ The 
presence of a second band gives an additional degree of 
freedom - in the sense of the previous paragraph - mak- 
ing the full many-particle problem very rich in new phe- 
nomena, yet also very complicated in general. We show 
how to effectively map the original problem to a much 
simpler spin system for specific fillings and parameters. 
This new model allows us to derive an analytical formula 
for the revival time of the interband oscillations which 
decay and revive due to weak inter-particle interactions. 



II. THE SYSTEM 

A. The many-body model 

We consider a two-band Bose-Hubbard model with 
an additional external force as obtained from a general 
many-body Hamiltonian under the assumption of a con- 
tact interaction and introduced in [3 [l9j]. We mea- 
sure all parameters in recoil energies E xec = h 2 k L /(2m), 



where kh is the wave vector of the laser creating the op- 
tical lattice and m the mass of the atoms. Setting h = 1 
throughout, the Hamiltonian reads [l8l.[l9| 
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The operator a; (at) annihilates (creates) a particle at 
site / of totally L sites in the lower band and bi (bj) in 
the upper band with the number operators nf — ajci^ 
n\ = b\b t The bands are separated by a bandgap A 
and have on-site energies ef — ±A/2 + IF. The hop- 
ping amplitudes between neighbouring sites in band a, b 
are denoted by t a , % > 0, and a repulsive interaction be- 
tween particles occupying the same site in band a (6) 
with a strength W a (Wb) has been included. The single- 
particle coupling of the bands is proportional to the ex- 
ternal Stark force F via CqF with a coupling constant 
Cq depending on the depth of the lattice Vo [TH . The 
bands are additionally coupled via the inter-particle in- 
teractions with a strength W x . Focusing on a realisation 
with a single optical lattice, the parameters fulfill gener- 
ally: A ^> t a ,tb, as well as t a ,tb « Wi and Co « —0.1. 
We take the external force F as a free parameter. Addi- 
tionally we assume that the interaction strength can be 
tuned, e.g., by the use of Feshbach resonances Q, and 
include a factor g to all interaction terms. For numeri- 
cal simulations, we change to the interaction-picture with 
respect to the external force HSl which removes the tilt 



lnf' b F and replaces a 
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e a i + i a i (and likewise 
for &J, . The Hamiltonian is then time-dependent with 
a periodicity of Tg = 2ir/F and allows to use periodic 
boundary conditions a^+i = <xi and = b\. 

To study the interband transport, we prepare the sys- 
tem in an initial state |^(0)), with a uniform distribution 
of particles in the lower band only and evolve it in time 
by the many-body Schrodinger equation. The quantity 
we study is the (normalised) number of particles in the 



2 



upper band 
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where N = + n i) i s the total number of parti- 

cles. We will refer to A/t>(£) as occupation of the upper 
band. For the range of parameters described above, this 
observable shows a superposition of many sinusoidal os- 
cillations with an amplitude of few per cent, even for 
strong forces [2l|. 



B. Weakly interacting system in resonance 

Despite the small interband oscillations described in 
the previous paragraph, a strong enhancement of the in- 
terband transport is possible for specific parameter val- 
ues. When the force-induced tilt of the lattice is such 
that a lower and upper band energy level become nearly 
degenerate (i.e. for A » mF, rn € N), the interband 
oscillations of Mb{t) come close to 100% indicating a res- 
onantly enhanced interband transport [2l[. We refer to 
these specific parameter values as resonant and will focus 
on this resonant behaviour in the following. 
To describe the non- interacting system H(g = 0) = Ho 
in resonance, it is useful to apply the following basis 
transformation involving Bessel functions of the first kind 
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with Xi — h/F, i = a,b. By using 
J2iez Jn-i{x)J n '-i{x') = J n -n>{x - x'), one finds 
that the transformation removes the hopping terms from 
the original Hamiltonian ([T]) but leads to a coupling 
between any sites from the two different bands weighted 
by Bessel functions 
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where Ax = x a + Xb and ef = ±A/2 + IF as above. The 
resonance condition, A « mF or equivalently e7 m er_ m) 
means that two levels from different bands become ener- 
getically degenerate. In this case, it is sufficient to keep 
only the direct coupling between these two sites leading 
to a sum of independent two-level systems 
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(5) 

This approximate description of the system in resonance 
corresponds to lowest order nearly degenerate perturba- 
tion theory and higher order corrections are easily calcu- 
lated, see e.g. [23|,[24j]. However, the lowest order approx- 
imation, eq. ([3]), gives already an accurate description of 



the non-interacting system in resonance [2lj, l23[. The 
resonant contribution to the non-interacting case is thus 
well-described as Mb(t) = sin 2 [CqF J m (Ax) t] with a pe- 
riod Xi cs = n/[CoFJ m (Ax)] ^> Tb much larger than the 
Bloch period. 

As demonstrated in [2l| , the inclusion of a weak inter- 
particle interaction leads to a dephasing of the resonant 
interband oscillations. The occupation of the upper band 
as a function of time exhibits a collapse and revival effect, 
with the time-scales for the collapse and revival inversely 
proportional to the interaction strength g [23[ . An exam- 
ple of such oscillations under a weak repulsive interaction 
is given by the solid line in fig. [2] below. In the weakly 
interacting regime under consideration here, one of the 
interaction terms in the full Hamiltonian eq. ([T]) is most 
important. We focus solely on repulsive interactions, for 
which the system tries to avoid double occupancy of sites 
in either bands. However, the system is always assumed 
to be at approximately integer filling and at the same 
time sites from different bands are nearly degenerate, 
such that it cannot avoid to have two particles occupying 
the same site in either band. Thus the dominant contri- 
bution comes from the term 2gW x ^ z n^n\ (see [2l[ for 
further details) . In the next section we derive an effective 
Hamiltonian that allows to study the effect of a weak in- 
teraction on the resonant interband oscillations in detail. 



III. RESULTS 

A. Effective spin model for system in resonance 

The description of the non-interacting system in reso- 
nance according to eq. §5$ contains already the seed for an 
effective model. The sum of many independent two-level 
systems can be viewed as a system of non-interacting 
spins. We only need to re-order the labeling of lattice 
sites such that the two levels being coupled have the same 
site-index and the coupling operator is then proportional 
to the Pauli matrix a x . The constant of proportionality 
is the coupling matrix element CqF J m {Ax) from eq. l[5]). 
This is simply a different way of writing the approximate 
Hamiltonian for the non-interacting system in resonance, 
eq. ([5]), and is schematically displayed for a resonance of 
order to = 2 in fig. [TJ To include the effect of the most 



important interaction term 2gW x nfrii into our effec- 
tive spin model, we insert the basis transformation from 
eq. ([3]) and obtain in the transformed basis 
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FIG. 1: Schematical view of the effect of the interband inter- 
action on the system in resonance of order m = 2. In reso- 
nance, the system is forced into a superposition of states from 
both bands (the sites forming a superposition are indicated 
by dashed ellipses). This happens on all lattice sites and the 
system cannot avoid an interaction of 2gW x Xlj n i n \- A new 
fictitious lattice labeling scheme setting the superpositions on 
one site is also indicated in the figure. 

Here we used the fact that only one of the many differ- 
ent combinations of Bessel functions gives a significant 
contribution [U [23[ . The reason is that the arguments 
of the Bessel functions x a and Xb are much smaller than 
unity for a realisation with a single optical lattice as dis- 
cussed here, and the dominant contribution is therefore 
given by the product of four zeroth order Bessel func- 
tions Jq (x a )J§ (xb)- We denote the interaction strength 
for this dominant process by 

U = 2gW x Jl{x a )Jl{x b ). (7) 

The introduction of the new ficticious lattice is now ef- 
fectively achieved by replacing I — y I + m for the sites of 
the upper band. It is important to note that the inter- 
band interaction was between atoms occupying the same 
site in different bands in the original lattice, i.e. oc n"n^, 
whereas in the new lattice it connects a particle at one 
lattice site in the lower band aJ a, = nf with a particle 
at a different site in the upper band Pi +m /3 l+m = rc£ , m 
(where we used the transformed basis ai,/3i). We fo- 
cus on unit filling N — L and, since the repulsive inter- 
action effectively suppresses higher occupation of lattice 
sites, we limit the occupation numbers of n? and nf to 
or I for our effective model. This allows us to replace 
n"'^ — > at'^ with the projectors on a spin-up or spin- 
down state aj l = (1; ± of ) /2. 

Collecting all arguments, the effective Hamiltonian (for 
a resonance of order m) is accordingly given by 

L 

H cS = (Vm<rf + Uajaf +m j (8) 
i=i 

where V m — CoFJ m (Ax). Here a\ denotes the Pauli 
matrices for a spin at site I. The first part is as in the 
non-interacting resonant system, which was also a sum 
of independent two-level systems. We only changed the 
ordering of the sites to bring degenerate levels close to- 
gether. The second part reflects the repulsion of two 
particles when sitting in different bands or different spin 
states respectively. Since we are using spin- 1/2 matri- 
ces in this effective description it can only be applied to 



the case of unit filling and the number of lattice sites 
is per definition identical with the number of spins. We 
expect it to be a good approximation for close-to- unit fill- 
ing (as is supported by our results below, see fig. 2]). An 
extension to higher fillings should be possible by using 
larger spins than spin 1/2, since this would allow fur- 
ther distinction of the type non-occupied, partly occu- 
pied, or highly-occupied, but is beyond the scope of the 
present article. The effective Hamiltonian ([5J is trans- 
lational invariant as our original model, such that one 
could use a reduction to subspaces of fixed total quasi- 
momentum similar to [l8|, HfJ, [23[. Please note, that 
the Hilbert space for the effective Hamiltonian has a di- 
mension dimH c ff = 2 L that is much smaller than the 
Hilbert space of the original bosonic problem, eq. ([lj, 
where dimH = (N + 2L - l)!/[JV!(2i - 1)!]. This is ad- 
vantageous for numerical computations since much larger 
system sizes become computable as compared to the orig- 
inal model. 

Let us discuss the effective model of eq. (|8j in more 
detail. The parameters in the effective spin model are 
chosen for the particular case of the system in resonance 
of order m. It includes only the resonant coupling be- 
tween the two sites and other non-resonant couplings are 
neglected. The effective model does thus not reproduce 
small scale oscillations which are found on top of the 
resonant oscillations within the full model (see [2l| for 
an example). However, these oscillations are only weakly 
influenced by the interparticle interaction and are not rel- 
evant for the collapse and revival effect we want to study. 
Another important aspect of the effective spin model con- 
cerns the choice of the interaction term a\ oy +m • Here we 
included only one of the four interaction terms from the 
original Hamiltonian, eq. (fT]), which has the strongest ef- 
fect on the resonant interband oscillations. The inclusion 
of the other terms is straightforward but not necessary 
in the present context. Furthermore we are limited to 
weak inter-particle interactions, U -C V m , which is, how- 
ever, not a limitation of the effective model but origi- 
nates from the physics of the original system: As soon 
as the inter-particle interaction becomes too strong, the 
resonant tunneling is washed out and the Rabi-like inter- 
band oscillations cede and eventually transform into an 
essentially structureless evolution of the band population 
defined in eq. © |23j]. 

To compare the effective model to the full problem, 
we computed the time-evolution of similar initial states 
in both models and show the resulting occupation of the 
upper band as a function of time with the pronounced col- 
lapse and revival effect in fig. @. The occupation of the 
upper band for the full model is given by Mb (t) , as defined 
by eq. ([2]), and has been computed by direct numerical 
integration of the time-dependent Schrodinger equation. 
In the effective spin model, a state with an atom occupy- 
ing the upper band is represented by a spin-up such that 
the corresponding observable for the spin model is given 
by Aff(t) — j;(ip(t)\ 0j IVK*)) an d the initial state is 
of the form LJ. ... 1) . Both observables are compared 
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FIG. 2: Occupation of the upper band in the full two-band 
Bose-Hubbard model (solid line; see [2ll ] for details) and num- 
ber of up-spins in the effective spin model (dashed line) for the 
same parameters. The spin model reproduces most features 
of the time signal, first and foremost the collapse and revival 
are in very good agreement between both models. Parame- 
ters correspond to Vo = 10: A = 7.77, t a = 0.005, t b = 0.121, 
Co = -0.114, W a = 0.040, W b = 0.027, and W x = 0.018; res- 
onance order m — 1, i.e. F — 7.9804; g = 0.1 and TV = 5 — L. 



in fig. [5] for a weakly interacting system of medium size. 
Overall, the effective spin model reproduces the occupa- 
tion of the upper band very well, especially when com- 
pared to the drastic simplification from the full two-band 
Bose-Hubbard model to the effective model of eq. ©. 
The good agreement is particularly surprising when tak- 
ing into account that a reduction of occupation numbers 
to or 1 is usually known as "hard-core bosons" [25|, [2(| 
and valid in the limit of strong interactions, whereas we 
are operating in exactly the opposite regime of U <C V m . 
Furthermore, the effective model reproduces the collapse 
and revival effect only when it is introduced from the 
transformed basis. Limiting the occupation number in 
the original a;,6;-basis by artificial constraints (such a 
truncation procedure wa s ap plied, e.g., in [27], HH) can- 
not reproduce the effect [231 ]. 

A great advantage of the effective spin model for the in- 
terband transport, eq. ©, is its exact solvability. Rewrit- 
ing the spin-up and -down operators in terms of Pauli ma- 
trices and applying a rotation of 7r/2 around the y-axis 
(which leads to a x — > a z and a z — > —cr x ), our effective 
Hamiltonian takes the following form 
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(V m af ~ \U ofoUm) + COnst - 



(9) 



This Hamiltonian is known (for m = 1) as the quantum 
Ising model in a transverse magnetic field [29[. It de- 
scribes coupled spins that tend to align in ^-direction but 
are subjected to the force of an applied magnetic field in 
x-direction. It can be solved exactly by subsequent appli- 
cation of a Jordan- Wigner transformation [30(, Fourier 
and Bogolyubov transformation. The final result in 



5(fc)(44-l/2). (10) 
The exact dispersion relation e(k) is given by 



k 



e(k) = 2V m ^l - ^cosfc + (^) 2 ps 2V m - \\J cos/c 

(11) 

and can be approximated for our weakly interacting sys- 
tem U «C V m as shown. Equation (fl0|) is the exact 
solution to our effective spin model. The elementary 
excitations of the system are non-interacting fermions 
with a dispersion relation that is approximately given 
by a cosine. These elementary excitations correspond to 
magnons, i.e., to delocalised spin-flips in the original spin 
basis. They read explicitly 

d k = cos{9 k /2)c k - i sin(6» fe /2)cl fc (12) 
sin k 

(13) 



tan 8 k 



cosk-Wm/U' 



where c k is the Fourier transform of c ; = a l e 11 ^ 11 "' 
with af = (erf ± iof )/2 and erf = 2c|q - 1. 



B. Revival time within the effective model 



The exact solution eq. (|T0|) of the effective Hamilto- 
nian allows, e.g., the computation of various correlation 
functions. But in the present context, we are interested 
in the time-evolution of particular initial states 



\E n ), 



(14) 



where \E n ) are the eigenstates to the effective model, 
eq. ([TO]) , and c n — (E n \ip(0)}. In general the overlaps 
c n between the given initial state and all eigenstates are 
needed for the time evolution. Instead of an analytical 
derivation of the overlaps on basis of the Jordan- Wigner 
transformation, we adopt a numerical approach here in 
order to decide which of the magnon states are relevant. 
In detail, it is sufficient to know which eigenstates have a 
significant contribution to the time-evolution to estimate 
the revival time of the resonant interband oscillations. 
Fig.[3]shows the coefficients c n for a time evolution of the 
initial state |1 . . . \) sorted by their eigenenergies. The 
eigenenergies appear in several bunches corresponding to 
eigenstates with a different number of magnon excita- 
tions ranching from to L magnons. Additionally, the 
three coefficients with the largest amplitude have been 
marked by squares in fig. [3] We find that the largest 
coefficients in the eigenbasis expansion are from the en- 
ergetically lowest eigenstates from the central bunches of 
the spectrum. To be more specific, we found numeri- 
cally that the largest coefficients always come from the 
subspaces with M — 1, M, and M + 1 magnons, where 
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FIG. 3: Shown are the coefficients c„ = (E n \ip(Q)) for an 
eigenbasis expansion of the initial state IV'(O)) = 1 4- ■•■ V) 
versus the corresponding eigenenergies (+). The three largest 
coefficients are marked by (□) and highlight eigenstates with 
M — 1, M, and M + 1 magnons. Parameters: V m = 1, L = 7, 
U = 0.25, and F — 4.6020. Note the logarithmic scale on the 
y-axis. 



M = L/2 for L even and M = (L - l)/2 for L odd, and 
are from the eigenstates with lowest energy within these 
subspaces. This important observation allows a simple 
estimate of the revival time as the time for a beating 
between oscillations with these three energies as frequen- 
cies. 

To find the lowest eigenenergy of state with M 
magnons, we use the fact that the energy of a many- 
body state with M magnons in the weakly interacting 
regime is according to eq. (jTTJ) given by 



M 

E M = ( 2Vm ~ \Ucos{k n )) , (15) 
i=i 



where kj t — 2irji/L and each ji can take a value between 
1, . . . , L. The energies for a given number of magnons 
M thus arise from different choices of the momenta kj t . 
The state with lowest energy in this cosine dispersion is 
obtained by using momenta that fill the empty cosine dis- 
persion from zero upwards. A many-magnon state with 
M = L/2 is reached when half of the possible states are 
filled and with Mil by adding or removing one magnon, 
respectively. This determines the momenta kj l to obtain 
a state with M magnons and minimal energy. We can 
now estimate the revival time from the difference between 
the energies of states with M— 1, M, and M+l magnons, 
i.e., we need Acj = (E M +i - E M ) - (E M - E M -i) = 
Em+i + Em-i — 2Em- Inserting explicitly that the en- 
ergy of M magnons is proportional to the sum of M 
cosine functions, eq. (|15[) . with different momenta filling 
the possible magnon states from below, we obtain the 



following frequency difference 

jj M+l M-l M 

Aw = - — ( ^ cos(%) + cos i k j) - 2 J2 cos ( k j)) 

i=i i=i j=i 

U 7T 

= -— (cosfc M+1 - cosk M ) m-—U, (16) 

where we expanded the cosine to lowest order around its 
zero. Using eq. (fT6|) and our expression for U, eq. (J7J, 
we find that the revival time as estimated by oscillations 
between the dominant frequencies is given by 

L 4tt 

The effective spin model predicts the revival time to be 
inversely proportional to the interaction strength and to 
a product of two Bessel functions from the basis transfor- 
mation. The parameters of the original full Hamiltonian 
eq. ((T|), like the hopping strengths, the gap between the 
two energy bands and the order of the resonance, en- 
ter via the arguments of Bessel functions x a ^ = t a j,/F, 
where the force has to be chosen according to the or- 
der of resonance F « A/to. These parameters and the 
revival time change when the depth of the optical lat- 
tice Vq is varied (see eq. @). Furthermore, the result 
eq. (|17|) from the effective spin model additionally pre- 
dicts a linear dependence of the revival time on the num- 
ber of lattice sites. In this way, our effective model adds 
an additional factor of L/2tt to our earlier result [2lj |. 
which has been obtained following the arguments that 
led to eq. ((TJ) above, with the assumption that the initial 
state is comparable to a coherent state and estimating 
the revival time by computing the effect of the dominant 
interaction term on a coherent state perturbatively [21| . 

To compare the result from eq. (|17p to numerical simu- 
lations of the full two-band Bose-Hubbard model, we use 
the size-dependent prefactor to rescale numerical results 
for different system sizes. The curves for different system 
sizes should coincide, as is verified in fig. 2J The revival 
times from full many-body models with Hilbert spaces 
ranging over three orders of magnitude fall onto one curve 
and demonstrate the validity of the effective spin model 
eq. ||5J) for the interband transport in the weakly inter- 
acting two-band Bose-Hubbard model. The remaining 
fluctuations with the system size originate in the approx- 
imation of eq. (|16l) and decrease with growing L. The 
additional offset arises from taking only the three largest 
coefficients for the derivation of the explicit expression 
in eq. (jTTJ) . This leads to an underestimation by about 
10% of the time for the maximum in the revived inter- 
band oscillations (which is our definition for the revival 
time). Inclusion of more than three coefficients should 
remove this systematic offset between the predicted and 
measured revival time. This has been tested by com- 
paring the oscillation dynamics of the original model ([1]) 
with the full spin model © [23J , and the good agreement 
is shown for an exemplary data point in fig. (cross, ex- 
tracted from the temporal evolution presented in fig. [5] 
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FIG. 4: We show the revival times from different numerical 
simulations of the full two-band Bose-Hubbard model (the 
corresponding system sizes 'AinL' and the dimension of the 
Hilbert spaces D are indicated in the legend). The numeri- 
cal values have been rescaled by the size-dependent prefactor 
2-k/L. The data points for different system sizes coincide, 
with small fluctuations due the approximation in eq. (|16|l . 
Our analytical prediction eq. (1 1 T|) (red dashed line) captures 
nicely the scaling with the system parameters, yet it shows a 
systematic offset arising from the three level approximation. 
The numerical result of the full spin model (cross marked by 
arrow), extracted from fig. [2] coincides with the full problem 
of eq. CEJ. 



model of eq. ([9]) for filling factors of the order one and 
for resonant coupling between the two energy bands. 
For weak inter-particle interactions the model is in good 
agreement with the full problem and allowed us to derive 
an analytical formula for the revivals of the resonant in- 
terband oscillations. Interesting future aspects to work 
on would be to include decay to higher energy states in 
the continuum part of the spectrum (e.g. by opening the 
model in a similar way as exercised in 



31] for a one- 



li 



band problem) and to extend the problem to atoms with 
internal structure. The internal degrees of freedom would 
become correlated with the external transport in "hori- 
zontal" - along the lattice - and "vertical" - between 
the bands - direction. Reductions of complex models are 
in general a necessary prerequisite in order to describe 
quantum systems with many degrees of freedom - pos- 
sibly of different kind and nature. So we hope that the 
spirit of our approach may inspire future research in this 
direction. 



above). We finally note, that eq. (fT7|) predicts a diver- 
gence of the revival time whenever the parameters of the 
system are chosen such that one of the Bessel functions 
in the denominator vanishes. This can be achieved, e.g., 
by tuning the energy gap A between the two bands, and 
the expected divergence of the revival time close to Bessel 
zeros was observed in numerical simulations [23[, giving 
room for a great deal of control of the resonant interband 
oscillations. 



IV. SUMMARY 

We have shown how to reduce the complexity of the 
original Hamiltonian of eq. ([!]) to the exactly solvable 
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